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Abstract 

Causal Dynamical Triangulations in four dimensions provide a background- 
independent definition of the sum over space-time geometries in nonperturbative 
quantum gravity. We show that the macroscopic four-dimensional world which 
emerges in the Euclidean sector of this theory is a bounce which satisfies a semi- 
classical equation. After integrating out all degrees of freedom except for a global 
scale factor, we obtain the ground state wave function of the universe as a function 
of this scale factor. 
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1. Introduction 

One important application of any theory of quantum gravity is a description of 
the quantum evolution of the very early universe. This is also the realm of quan- 
tum cosmology, which tries to capture the essence of the gravitational dynamics 
by quantizing only a finite number of degrees of freedom characterizing the uni- 
verse as a whole. The path integral formulation of quantum cosmology came 
to prominence with the work of the Cambridge group and others on Euclidean 
quantum gravity [T|, and in particular that of Hartle and Hawking [2J. Central 
in this and related approaches is the construction of a "wave function of the uni- 
verse" , either as a solution of the Wheeler-DeWitt equation or a propagator for 
the theory (see, for example, |I3 EJ lEJ EH E|). In attempting to do this, a variety 
of technical and conceptual issues has to be addressed, including the choice of 
boundary conditions for the wave function, the unboundedness of the gravita- 
tional action and ensuing divergence of the Euclidean cosmological path integral, 
the appropriateness of the minisuperspace and/or semiclassical approximations, 
and the physical interpretation of the construction (see [H] for a recent concise 
review) . 

One could hope that a nonperturbative path integral formulation which does 
not impose any a priori symmetry restrictions on the geometry of the universe 
would help resolve some of these issues. "Causal Dynamical Triangulations" 
provide exactly such a background-independent, nonperturbative definition of 
quantum gravity, in which the sum over all space-time geometries is constructively 
defined and the causal (Lorentzian) structure of space-time plays a crucial role 
[3 GH El El It can be viewed as a realization of an idea of Teitelboim's, 
who argued that in a (continuum) proper-time formulation of the Lorentzian 
gravitational path integral one should integrate over positive lapse functions only, 
thereby building a notion of causality into the quantum dynamics |14j . 

In the context of quantum cosmology it has been argued [UJ El that a tun- 
neling wave function a la Vilenkin, where a universe tunnels from a vanishing to 
an extended three-geometry, is a special case of Teitelboim's causal propagator 
between two three-geometries. In the present work, where the wave function of 
the universe will be constructed from first principles, we will indeed observe a 
similar phenomenon, although our interpretation in the end will be somewhat 
different. 

In [T3| we reported that the approach of Causal Dynamical Triangulations, de- 
spite its background independence, generates a four-dimensional universe around 
which (small) quantum fluctuations take place. The purpose of this letter is to 
identify the effective action which determines the shape of this macroscopic 4d 
world. Rather surprisingly we find that the effective action which describes the 
infrared, long-distance part of the universe is closely related to a simple minisuper- 
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space action frequently considered in quantum cosmology. The only differences 
in our full quantum treatment are that (i) the unboundedness problem of the 
conformal mode in the Euclidean sector is cured, and (ii) the ultraviolet, short- 
distance part of the effective action is such that the solution to the Euclidean 
action describes a bounce from a universe of no spatial extension to one of finite 
spatial size. While resembling Vilenkin's picture of a "universe from nothing" [3], 
the interpretation in the present context is rather in terms of the ground state 
wave function of the universe with everything but the scale factor integrated out. 
We describe how one can determine this wave function from first principles. 

It should be emphasized that unlike in standard minisuperspace models for 
cosmology, we do not assume homogeneity or isotropy, nor do we impose any other 
a priori symmetry conditions on the gravitational degrees of freedom. We perform 
the full path integral and determine the effective Lagrangian which describes the 
dynamics of the global scale factor, as well as the ground state wave function of 
the universe as a function of this scale parameter. 

The rest of this article is organized as follows: in the next section we recall 
some salient features of the Causal Dynamical Triangulations approach, includ- 
ing the set-up of the numerical simulations and recent numerical results in four 
dimensions. We then demonstrate in Sec. 3 that the numerical data are perfectly 
described by a simple minisuperspace action. Sec. 4 outlines how this result re- 
lates to the ground state wave function of the universe, and the final Sec. 5 is 
devoted to a discussion. 

2. Observing the bounce 

The idea to construct a quantum theory of gravity by using Causal Dynamical 
Triangulations was motivated by the desire to formulate a quantum gravity theory 
with the correct Lorentzian signature and causal properties |14j . and to have a 
path integral formulation which may be closely related to attempts to quantize 
the theory canonically. For the purposes of this letter, we will only summarize the 
main properties of this approach; more details on the rationale and techniques 
can be found elsewhere jHJ EH HU H2j . 

We insist that only causally well-behaved geometries appear in the path in- 
tegral, which is regularized by summing over a particular class of triangulated, 
piecewise flat (i.e. piecewise Minkowskian) geometries. All causal simplicial space- 
times contributing to the path integral are foliated by a version of "proper time" 
t, and each geometry can be obtained by gluing together four-simplices in a way 
that respects this foliation. Each four-simplex has time-like links of negative 
length- squared — and space-like links of positive length-squared a^, with all of 
the latter located in spatial slices of constant (integer in lattice units) proper time 
t. These slices consist of purely space-like tetrahedra, forming a three-dimensional 
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piecewise flat manifold, whose topology we choose for simplicity to be that of a 
three-sphere S 3 . A necessary condition for obtaining a well-defined continuum 
limit from this regularized setting is that the lattice spacing a oc a t oc a s goes to 
zero while the number N4 of four-simplices goes to infinity in such a way that 
the continuum four- volume V4 := a 4 A^ stays fixed. Let us emphasize that the 
parameter a therefore does not play the role of a fundamental discrete length. 
A further property of our explicit construction is that each configuration can be 
rotated to Euclidean signature, a necessary prerequisite for discussing the con- 
vergence properties of the sum over geometries, as well as for using Monte Carlo 
techniques. 

The partition function for quantum gravity is 

Z(A,G) = J V[g] S[g] = ^J #xy/\6&g{ (R - 2A), (1) 

where S[g] is the Einstein-Hilbert action including a cosmo logical-constant term 
A, and G the gravitational constant 1 . Using our simplicial regularization this 
becomes 

Z(A,GW = £-L e ^l, ( 2 ) 
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where the integration over Lorentzian geometries is replaced by a sum over causal 
triangulations, and Ct is a symmetry factor of the triangulation T, the order of its 
automorphism group. The action S[T] is Regge's version of the Einstein-Hilbert 
action, appropriate for piecewise linear geometries (see J2| for details). In the 
remainder of this article we will for simplicity use a continuum notation, but 
it should be understood that whenever computer simulations are mentioned the 
implementation of the path integral is in terms of piecewise flat geometries. 

The substitution — a\ — > turns all time-like into space-like edges and rotates 
all configurations g to Euclidean space-times gg, replacing at the same time 
iS[g] — ► —Se[9e}i where Se denotes the Euclidean Einstein-Hilbert action. In 
computer simulations it is often convenient to work with universes of constant 
four- volume V4. The Euclidean counterpart of the partition function (JIJ can be 
decomposed as 

poo 

Z E (A,G)= / dV 4 e-^ AV *Z E (V 4l G), (3) 
Jo 

where the partition function Ze(Vi,G) for fixed four-volume is defined as 

Z E (V A ,G)= Jv[g]e- §E ^8(Jd i xy^tg'-V A ), S E [g] = ~ J d i xy^tg' R. 

(4) 



1 We ignore a numerical constant multiplying G in QJ. 
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Whenever V4 is kept fixed we will use Z E (V±, G) as our partition function. It is 
related to Z#(A, G) by the Laplace transformation (jSJ). 

The specific (Euclidean) partition function we will consider is the so-called 
quantum-gravitational proper-time propagator defined by 



where the integration is over all four-dimensional (Euclidean) geometries g& of 
topology S 3 x [0, 1], each with proper time running from to t, and with spatial 
boundary geometries #3(0) and g 3 (t) at proper times and t. 




Figure 1: Monte Carlo snapshot of a "typical universe" of discrete volume 91.100 
four-simplices and total time extent (vertical direction) t = 40. The circumference at 
integer proper time s is proportional to the spatial three- volume V^s). The surface 
represents an interpolation between adjacent spatial volumes, without capturing the 
actual 4d connectivity between neighbouring spatial slices. 

While it may be difficult to find an explicit analytic expression for the full 
propagator (jSJ) of the four-dimensional theory, Monte Carlo simulations are read- 
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ily available, using standard techniques from Euclidean dynamically triangulated 
quantum gravity ^E]. For convenience of the computer simulations, we keep the 
total four-volume V4 of space-time fixed and also often use periodic rather than 
fixed boundary conditions, i.e. sum over space-times with topology S 3 x S 1 rather 
than S 3 x [0, 1]. This periodicity does not affect the results reported below, as is 
illustrated by FigJTJ This shows the typical "shape" (spatial three- volume V 3 (s) 
as a function of proper time s) of a space-time configuration generated by the 
computer. 2 For given space-time volume V4, as long as t is chosen sufficiently 
large, the configuration will develop a thin stalk like the one shown in Fig.l. It 
will then not matter for the analysis of the large-scale geometry whether or not 
time is periodically identified. 

A convenient "observable" is the spatial volume- volume correlator defined by 

CV 4 (A) = (V 3 (0)V 3 (A)) V4 = ^l ds (Ms)V 3 (s + A)) V4 , (6) 

where we have identified the spatial boundary geometries at times and t. We 
have measured this correlator for various four-volumes V4 and plotted the result 
as a function of the scaled variable x = A/V^ d , after normalizing the correlator 
for a given V4 to have an integral equal to 1 (in lattice units). The variable d is 
then determined from the condition that the overlap between the correlators of 
different V4 be maximal. 3 Fig |2] illustrates the almost perfect overlap obtained for 
d = 4. We take this as strong evidence that our "macroscopic" space-times (we 
are using up to 360.000 four-simplices) are genuinely four-dimensional. 4 As em- 
phasized in ^3] > this is a highly non-trivial result and does not follow from the fact 
that the individual building blocks at the cut-off scale a are four-dimensional. Ad- 
ditional evidence for a macroscopically four- dimensional universe was presented 
in jinj, and an extended analysis will appear in due course [T%] . 

Our goal here will be to understand the precise analytical form of the volume- 
volume correlator CV 4 (A). To this end, let us consider the distribution of differ- 
ences in the spatial volumes V3 of successive spatial slices at proper times s and 
s + 5, where 5 is infinitesimal, i.e. 5 = 1 in lattice proper time units. We have 
measured the probability distribution Py 3 (z) of the variable 

g= r,(. + f)-tW V3 = y 3(s) + V3{B + S y (7) 

^3 

2 "Typical" in this context is used to characterize the geometric characteristics that will be 
shared with probability 1 by a randomly chosen member of the ensemble of four-geometries. 

3 We are employing standard finite size scaling methods from the theory of critical phenomena 
in statistical mechanics, see, for example |17j . 

4 Similar data, but with a smaller maximal four- volume, were already published in |Ki|. 
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Figure 2: Measurement of spatial volume- volume correlator for space-times with 
22.250, 45.500, 91.000, 181.000 and 362.000 four-simplices, plotted as function of the 
scaled variable x = A/V^ 4 . 

for different values of V3. As shown in Fig. |3] they fall on a common curve. 5 
Furthermore, the distribution Pv 3 (z) is fitted very well by a Gaussian e~ cz , with 
a constant c independent of V3. From estimating the entropy of spatial geome- 
tries, that is, the number of such configurations, one would expect corrections of 
the form V3*, with < a < 1, to the exponent cz 2 in the distribution Py 3 (z). 
Unfortunately it is impossible to measure these corrections directly in a reliable 
way. We therefore make a general ansatz for the probability distribution for large 
V 3 (s) as 

exp 

where < a < 1, and c\ and C2 are positive constants. 

We are thus by "observation" led to the following effective action for large 
three- volume V 3 {s): 

5 Again we have applied finite size scaling techniques, starting out with an arbitrary power 
Vjj Q in the denominator in l(7)l. and have determined a = 1/2 from the principle of maximal 
overlap of the distributions for various V^'s. 
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Figure 3: Distribution Py 3 (z) of volume differences of adjacent spatial slices, for three- 
volumes V 3 = 10.000, 20.000, 40.000 and 80.000 tetrahedra. 

where A is a Lagrange multiplier to be determined such that 

/ dsV 3 (s) = V 4 . (10) 
Jo 

From general scaling of the above action it is clear that the only chance to obtain 
the observed scaling law, expressed in terms of the variable t/V^ 4 , is by setting 
a = 1/3. In addition, to reproduce the observed stalk for large times t the 

1/3 

function V 3 has to be replaced by a function of V 3 whose derivative at goes 
like V 3 U , v > 0, for reasons that will become clear in Sec. 3 below. A simple 
modification, which keeps the large- V3 behaviour intact, is given by 

^ 3 1/3 ^(l + ^) 1/3 -l, (11) 

but the detailed form is not important. If we now introduce the (non-negative) 
scale factor a(s) by 

V 3 (s) = a 3 (s), (12) 
we can (after suitable rescaling of s and a(s)) write the effective action as 

se " = ^l ds ( a(s) (^r) 2 + a(s) " Afl3(s) ) ' (13) 

with the understanding that the linear term should be replaced using (|12|) and 
(fTTj) for small a(s). We emphasize again that we have been led to ()13j) entirely 
by "observation" and that one can view the small-a(s) behaviour implied by (jllj) 
as a result of quantum fluctuations. 
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3. Minisuperspace 

Let us now consider the simplest minisuperspace model for a closed universe 
in quantum cosmology, as for instance used by Hartle and Hawking in their 
semiclassical evaluation of the wave function of the universe j2J- In Euclidean 
signature and proper-time coordinates, the metrics are of the form 

ds 2 = dt 2 + a 2 {t)dtt 2 3 , (14) 

where the scale factor a(t) is the only dynamical variable and dOg denotes the 
metric on the three-sphere. The corresponding Einstein-Hilbert action is 

Seff = ^jdt (~a(t) 2 - a (t) + Aa 3 (t) j . (15) 

If no four- volume constraint is imposed, A is the cosmological constant. If the 
four-volume is fixed to V4, such that the discussion parallels the computer simu- 
lations reported above, A should be viewed as a Lagrange multiplier enforcing a 
given size of the universe. In the latter case we obtain the same effective action 
as in (fT3|) up to an overall sign, due to the infamous conformal divergence of the 
classical Einstein action. Let us for the moment ignore this overall minus sign and 
compare the two potentials relevant for the calculation of semiclassical Euclidean 
solutions associated with the actions (fT5|) and (fTB"|) . The "potential" 6 is 

V( a ) = -a + \a 3 , (16) 

and is shown in Fig. without and with small-a modification, for the stan- 
dard minisuperspace model and our effective model, respectively. The quantum- 
induced difference for small a is important since the action (|13|) allows for a classi- 
cally stable solution a(t) = which explains the "stalk" observed in the computer 
simulations. Moreover, it is appropriate to speak of a Euclidean "bounce" be- 
cause a = is a local maximum. If one therefore naively turns the potential 
upside down when rotating back to Lorentzian signature, the metastable state 
a(t) = can tunnel to a state where a(t) ~ V^ 4 , with a probability amplitude 
per unit time which is (the exponential of) the Euclidean action. We will discuss 
this further in the next section. 

In order to understand how well the semiclassical action ()13|) can reproduce 
the Monte Carlo data, that is, the correlator CV 4 (A) of FigHJ we have solved for 
the semiclassical bounce using (fT3|) . and presented the result as the continuous 

6 To obtain a standard potential - without changing "time" - one should first transform to a 
variable x — for which the kinetic term in the actions assumes the standard quadratic form. 
It is the resulting potential V(x) = — x 2 / 3 + \x 2 which in the case of (|13fl should be modified 
for small x such that V"'(0) = 0. 
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Figure 4: The potential V(a) of (|16|) underlying the standard minisuperspace dynamics 
(left) and the analogous potential in the effective action obtained from the full quantum 
gravity model, with small-a modification due to quantum fluctuations (right). 



black curve in FigJ21 7 . The agreement with the real data generated by the Monte 
Carlo simulations is clearly perfect. 

4. The wave function of the universe 

The picture emerging from the above for the effective dynamics of the scale factor 
resembles that of a universe created by tunneling from nothing (see, for example, 
IH El), although the presence of a preferred notion of time makes our situ- 
ation closer to conventional quantum mechanics. In the set-up analyzed here, 
there is apparently a state of vanishing spatial extension which can "tunnel" to a 
universe of finite linear extension of order a ~ V^ 4 . Adopting such a tunneling 
interpretation, the action of the bounce is 



and the associated probability per unit proper time for the tunneling given by 



However, recall that this picture arose from a situation where for computer- 
technical reasons we imposed a constraint on the four-volume. Since our formu- 

7 More precisely, we solved the classical equation of motion corresponding to the potential 
shown in Fig0]on the right, with an energy slightly below zero (the closer to zero the longer 
the stalk), and used this solution to create an artificial distribution of three- volumes Va(s) 
analogous to the one generated by Monte Carlo simulation from first principles. We then 
treated this artificial distribution precisely as if it had come from real Monte Carlo data. 




(18) 
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lation possesses a well-defined Hamiltonian 8 which is bounded below, the ground 
state wave function can be chosen real and positive. In view of this, calling (|18|) a 
tunneling probability is misleading, since it would imply an oscillating behaviour 
for a ^> V^ 1//4 . The correct interpretation of ()18|) is rather that of the square of 
the ground state wave function for a ~ V^ 4 . 

To illustrate what we have in mind, let us consider a quantum-mechanical 
system with Hamiltonian H = p 2 /2 + V(x), where the minimum of the potential 
V is at x = 0. The ground state wave function has the path integral representation 

/•x(0)=a 

*o(o) ~ / T>x(t) e~ SE[x{t)] , (19) 



J x(— oo)=0 

where «9e [#(£)] is the classical Euclidean action 

1 



S E [x(t)} = f dt 

J — oo 



-x 2 - (-V(x)) 



(20) 



If there is a classical solution x c i(t) which extremizes the Euclidean action (J20)) 
with boundary conditions x c /(— oo) = and x c ;(0) = a, a semiclassical calculation 
of ^o(a) involves a saddle point expansion around that solution and the leading 
exponential of the wave function is 



%{a)~e- SE ^W, S E [x d (t)}= / dx^W{x). (21) 



As an example, for the harmonic oscillator we have x c i(t) = ae"' and (J2T|) is 
exact. For a general potential the semiclassical approximation will of course not 
be exact. Nevertheless, we have presented strong evidence that in the case of 
quantum gravity, integrating over all degrees of freedom except the three-volume 
and defining a = V^ 3 , the semiclassical approximation is excellent. If we assume 
that it is equally good in the absence of the four- volume constraint, we compute 
in a straightforward way from (|21|) that 

* (a) ~ e -ro((i+A» 2 ) s/a -i) > (22) 

with c a constant of order one. Note that A in f!22D is the real cosmological 
constant and no longer a Lagrange multiplier. We have thus calculated the wave 
function of the universe from first principles up to prefactors and corrections to 
the semiclassical approximation. 

8 In the framework of Causal Dynamical Triangulations one has a transfer matrix between 
consecutive proper-time slices whose logarithm gives in principle the full quantum Hamiltonian, 
see ^2] for details. 
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It is important to understand that the wave function ^o(a) can be addressed 
via computer simulations using a decomposition analogous to (J3J), namely, 




where the functional integration in (|24J1 is over four-geometries with V^— oo) = 
and V^t — O) = a 3 . The computer simulations reported here were done for the 
special case V^{t = 0) = 0. 9 Using (}23|) and (J24|) one can now check whether the 
semiclassical approximation reported here is valid for all values of a and V4. If 
so, one will be led to (|2*2*|) . 

5. Discussion 

Causal Dynamical Triangulations constitute a framework for defining quantum 
gravity nonperturbatively as the continuum limit of a well-defined regularized 
sum over geometries. We reported recently on the outcome of the first Monte 
Carlo simulations in four dimensions l3j. Very encouragingly, we observed the 
dynamical generation of a macroscopic four-dimensional (Euclidean) world, with 
small quantum fluctuations superimposed. In this letter we showed that the scale 
factor characterizing the macroscopic shape of this ground state of geometry is 
well described by an effective action similar to that of the simplest minisuperspace 
model used in quantum cosmology. However, in our case such a result has for the 
first time - we believe - been derived from first principles. 

The negative sign of the kinetic term in the standard minisuperspace action 
(I15|) reflects the well-known unboundedness of the conformal mode in the Eu- 
clidean Einstein-Hilbert action. Amazingly, after integrating out all variables 
except the scale factor (which is simply the global conformal mode), we obtain a 
positive kinetic term in our effective action (|13|). This is consistent with a contin- 
uum formulation of the gravitational path integral in proper-time gauge, where 
a strong argument was made for the nonperturbative cancellation of the confor- 
mal divergence by a measure factor coming from a Faddeev-Popov determinant 
[19J. Similar ideas were pursued in [201 ; although in that case matter fields were 
needed to change the sign of the conformal mode term. The phenomenon of sign 
change of the conformal kinetic term is also familiar from the 2d Euclidean quan- 
tum theory where, due to the conformal anomaly, integrating out unitary matter 
with central charge < c < 1 yields an effective term — c(d<fr) 2 for the conformal 

9 Since the spatial volume of the universe is not a monotonic function of V4, the integration 
range in i|24l) in this case should be split into two intervals [0, a max ] and [a max ,0]. 
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factor in the action, making it unbounded below. Again it is a Faddeev-Popov 
determinant, arising from the requirement to integrate not over metrics, but only 
geometries, which adds a 26 (90) 2 and ensures that the combined kinetic term is 
positive. 

A number of open issues remain to be addressed, including the details of 
the renormalization mechanism. Here Causal Dynamical Triangulations gives 
us the possibility to study Weinberg's scenario of "asymptotic safety" [21] in 
the context of an explicit quantum-gravitational model. As indicated in earlier 
work on Causal Dynamical Triangulations in space-time dimension three, the 
renormalization may be non-standard [22], which in a way would be welcome. 
This is also supported by the present computer simulations, in the sense that no 
fine-tuning of the bare gravitational coupling constant seems to be necessary to 
reach the continuum limit. Details of this will be discussed elsewhere |18j . 

A most interesting question is of course how the above semiclassical cosmolog- 
ical picture is changed by the inclusion of matter fields. We have now the chance 
to investigate a number of possible scenarios suggested in quantum cosmology 
from first principles. 
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